Quantized spin wave modes in magnetic tunnel junction nanopillars 
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We present an experimental and theoretical study of the magnetic field dependence of the mode frequency of 
thermally excited spin waves in rectangular shaped nanopillars of lateral sizes 60 x 100, 75 x 150, and 105 x 
190 nm 2 , patterned from MgO-based magnetic tunnel junctions. The spin wave frequencies were measured 
using spectrally resolved electrical noise measurements. In all spectra, several independent quantized spin wave 
modes have been observed and could be identified as eigenexcitations of the free layer and of the synthetic 
antiferromagnet of the junction. Using a theoretical approach based on the diagonalization of the dynamical 
matrix of a system of three coupled, spatially confined magnetic layers, we have modeled the spectra for the 
smallest pillar and have extracted its material parameters. The magnetization and exchange stiffness constant 
of the CoFeB free layer are thereby found to be substantially reduced compared to the corresponding thin film 
values. Moreover, we could infer that the pinning of the magnetization at the lateral boundaries must be weak. 
Finally, the interlayer dipolar coupling between the free layer and the synthetic antiferromagnet causes mode 
anticrossings with gap openings up to 2 GHz. At low fields and in the larger pillars, there is clear evidence 
for strong non-uniformities of the layer magnetizations. In particular, at zero field the lowest mode is not the 
fundamental mode, but a mode most likely localized near the layer edges. 

PACS numbers: 75.75.+a, 75.30.Ds, 85.75.-d, 84.40.-x 



I. INTRODUCTION 

In the last few years, magnetic tunnel junction (MTJ) 
nanopillars have received tremendous attention due to their 
promising potential for applications in spin-transfer-switched 
Magnetic Random Access Memory or as spin-torque oscil- 
lators for microwave generation^ With GHz frequencies 
the operation speed of these devices happens to lie in the 
same frequency range as the dynamic eigenexcitations of 
the underlying nanoelements (thermally excited spin waves), 
which may therefore manifest themselves as unwanted noise 
sources. However, as eigenexcitations, thermal spin waves 
also constitute an excellent probe for the intrinsic magnetic 
properties of the nanopillars. The experimental detection of 
spin waves in MTJ nanopillar devices and the understanding 
of their nature is therefore of great interest both for funda- 
mental and technological reasons. 

Spin waves in confined structures have been studied 
extensively in single-layer dots with thicknesses between 40 
and 15 nm and typical lateral dimensions from 3 /im down to 
200 nm. 4 " 8 In these systems, two types of spin wave modes 
have been identified: quantized volume modes located around 
the center of the element where the internal field is basically 
homogeneous, and spin wave well or end modes localized 
near the element edges in the inhomogeneity region of the 
internal field. The above elements are characterized by their 
thickness being significantly larger than the exchange length 
of the layer material (typically 5 nm). In structures with this 
property the dominating interaction is the magneto-static 
dipolar interaction, 9 which causes the inhomogeneity of 
the internal field, thus determining the character and spatial 



profile of the modes. 

The eigenexcitations of a multi-layer dot differ in general 
significantly from those of an ensemble of isolated magnetic 
dots due to the interlayer interactions between the magnetic 
layers in the stack: mutual dipolar coupling and - for suf- 
ficiently thin metallic spacer layers - interlayer exchange 
coupling^ 

Eigenexcitations of nanopillar structures have been the sub- 
ject of very few studies so far. Thermal spin waves have been 
investigated systematically only in pseudo-spin-valves 1 ^— 
of circular and elliptical shape (smallest dimension 200 nm) 
consisting of two magnetic layers of 10 nm thickness 
separated by a 10 nm thick spacer layer, i.e. again layer 
thicknesses were much larger than the exchange length. 
Consequently, the profiles of the modes in each of the two 
pillar layers showed great resemblance 712 with the mode 
profiles in the corresponding isolated dots. In a symmetric 
spin-valve stack, 11 the main impact of the mutual dipolar 
coupling between the layers was found to be a fixed phase 
relation between the modes in the two layers for high applied 
field, and hybridization effects at low field. 

Common MTJ nanopillars differ qualitatively from the 
pseudo-spin-valves in three fundamental points: Firstly, 
with a free layer and a synthetic antiferromagnet (SAF) they 
consist of three magnetic layers; secondly, with 2 — 4 nm the 
layer thicknesses are now smaller than the exchange length, 14 
such that, for sufficiently small lateral dimensions, the spin 
dynamics in each layer is dominated by the exchange inter- 
action; thirdly, the interlayer interaction of the three layers is 
highly asymmetric: the two SAF layers are strongly coupled 
by interlayer exchange and - more weakly - mutual dipolar 
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coupling, one of them (the pinned layer) being additionally 
subject to the strong exchange bias field; the free layer 
interacts with the SAF via the comparatively weak mutual 
dipolar coupling only. Consequently, the eigenexcitations of 
an MTJ nanopillar are expected to be much more complex 
than those of the pseudo-spin-valves. 

In this paper, we investigate the magnetic field dependence 
of the mode frequency of thermally excited spin waves in rect- 
angular shaped MgO-based MTJ nanopillars of different lat- 
eral sizes. In section II, we will describe the basic magnetic 
properties of the devices and the experimental techniques used 
to acquire the spin wave spectra. The features of the measured 
spectra in dependence of the pillar size and the direction of the 
applied field are described in the following section III. In sec- 
tion IV, we will point out short-comes of the macrospin model 
when applied to our samples, as a consequence of which we 
will introduce in section V a model of quantized spin wave 
modes in nanopillars consisting of three magnetic layers. In 
section VI, we will use this model to extract the material pa- 
rameters of the pillar, which will finally be discussed in sec- 
tion VII along with the limitations of our model. 



II. SAMPLES AND EXPERIMENTAL TECHNIQUES 

A. Samples and basic device properties 

The fabrication and basic properties of our samples 
are described in Ref. they are rectangular shaped 

nanopillars, all patterned from the same MTJ stack of com- 
position Co6oFe2oB2o (3 nm, free layer)/ Mg(1.3)[nat. ox.]/ 
Co6oFe2oB2o(2, reference layer)/ Ru(0.8)/Co7oFe3o(2, 
pinned layer)/PtMn(20), deposited by Singulus Technologies 
AG. The three layers following the MgO tunnel barrier 
compose the synthetic antiferromagnet (SAF). The pillars 
were designed in three lateral sizes: 60 x 100, 75 x 150, 
and 105 x 190 nm 2 , which will be referred to as small (S), 
medium (M), and large (L) size, respectively. Note that 
unlike in Ref. [l5| the given dimensions are not the nominal 
values, but mean values measured on the exposed e-beam 
resist with a device-to-device deviation of ±10 nm. In order 
to obtain electrically contactable devices the nanopillars were 
inserted in series between coplanar leads, following design 
rules ensuring high bandwidth^ 

The devices have a resistance area product of 16 O /im 2 
and typically 80% tunnel magneto-resistance ratio. Their hys- 
teretic properties are consistent with the uniaxial anisotropy 
expected from the rectangular pillar shape where the long 
edge of the rectangle, oriented along the exchange pinning 
direction of the PtMn antiferromagnet, is the easy axis (EA), 
and the short edge the hard axis (HA) of the magnetization. 
Panels (b) and (d) of Fig. Q] show as a reference EA and 
HA hysteresis loops of a nanopillar of size S calculated 
in macrospin approximation using as material parameters 
literature bulk values (see figure caption). In comparison, the 
experimental EA and HA loops for devices of pillar size S, 



M, and L are depicted in panels (b), (d), and (f) of Figs. 12131 
respectively. 

At negative EA applied fields the pillars are in the parallel 
(P) state, at positive fields in the antiparallel (AP) state; 
spin-flop (SF) transition of the SAF occurs typically at EA 
fields around +170 mT. Room temperature coercivity is 
25 — 35 mT for devices of size S and M, and 20 — 25 mT 
for size L. From astroid measurements 15 mean anisotropy 
fields of 37 mT, 46 mT, and 38 mT for pillar sizes S, M, and 
L, respectively, have been determined. The EA loops of all 
devices are off-centered towards negative fields, indicating 
non-negligible antiparallel coupling of the free layer magne- 
tization and the SAF. With increasing pillar size this coupling 
is decreasing: while for size S the shift is 5 — 11 mT, it is only 
3 — 7 mT for size M, and 1 — 5 mT for size L. 
The bell shape of the HA hysteresis loops (Fig. [3]) is consis- 
tent with the antiparallel coupling observed on the EA: At 
zero HA applied field the devices are always in the AP state. 
With increasing (absolute value of the) field the resistance 
decreases continuously from the maximum resistance of 
the AP state down to almost parallel remanence, as the 
magnetizations of the free layer and finally the two SAF 
layers progressively tilt towards the applied field. Ascending 
and descending field branch of the HA loops are for most 
devices identical. In contrast, in the absence of coupling 
between free layer and SAF, the pillar relaxes with equal 
probability into P state or AP state when the HA field is 
switched off, resulting in two branches of the hysteresis loop 
(see Fig. |TJd)). Finally, note that the sharp bends in the 
resistance curve at about ±60 mT in the measured HA loop 
for pillar size L, and also present in the calculated HA loop, 
become more and more rounded for decreasing pillar size, i.e. 
for increasing antiparallel coupling. 

We have used the intrinsic symmetry of the HA loops to 
align the external field with the symmetry axes of the rectan- 
gle by choosing the field direction such that the loops showed 
highest possible symmetry. For some devices the loops were, 
though symmetric at high fields, noticeably asymmetric at low 
fields; in these cases we have cross-checked the alignment 
with the symmetry of the corresponding HA spectra. The mis- 
alignment of the field should therefore not exceed 2°. 



B. Set-up and experimental methods 

To obtain their spin wave spectra the devices were inserted 
into a high bandwidth circuit similar to that in Ref. 
and their voltage noise power spectrum density (PSD) was 
measured for moderate dc bias currents as a function of the 
applied magnetic field. The noise spectrum at each field 
step was obtained by subtracting from the spectrum for 
non-zero bias current a zero-current reference spectrum in 
order to eliminate noise of non-magnetic origin. In panels 
(a), (c), and (c) of Figs. 1213 1 examples of 2D density plots 
of the PSD versus the magnetic field are shown, where the 
dark regions correspond to maxima in the PSD and therefore 
to eigenexcitations of the magnetic system, 18 i.e. spin wave 
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modes. The spectra are displayed in a contrast scaling with 
the logarithm of the noise in excess to the noise at zero 
current. As the difference in amplitude between the most 
intense and the weakest modes is even on a logarithmic scale 
still large, in all figures the gray-scale of the PSD has been 
modulated, and black dots have been superimposed to better 
evidence the weaker modes. 

Bias currents used to measure the spectra were chosen as low 
as possible in order not to affect the mode frequencies, but 
still high enough to obtain sufficient signal-to-noise ratio. 
Devices of pillar size S were therefore mostly measured at 
±0.1 mA, those of size M at ±0.2 mA, and those of size L at 
±0.3 mA. The differences in amplitude for opposite current 
polarity were hardly noticeable on the logarithmic scale, and 
the maximum difference in the mode frequencies, observed 
for 0.3 mA, was 0.15 GHz. These observations are in 
agreement with previous works 19 on similar samples, where 
a spin-torque threshold current of 1 .6 mA for size-L devices 
has been determined. On the frequency scales considered 
in this paper the bias current dependence of the spectra can 
therefore be neglected. 

Note that our measurement technique allows to detect the 
spin wave modes of individual pillars, in contrast to the Bril- 
louin Light Scattering technique used in Refs. l4lfTl1fT3l or the 
frequency-domain coplanar waveguide technique of Ref. d, 
where the measured spectra were an average over a large num- 
ber of devices. Moreover, we do not need optical or any other 
direct access to the magnetic layers of the pillar, but can mea- 
sure them in their natural working environment, i.e. as part 
of a stack used in (actual or potential) functional devices, and 
subject to electrical currents. Finally, since we measure the 
voltage noise, which is in one-to-one correspondence with the 
magneto-resistance (MR) noise, we are equally sensitive to 
spin waves in the free layer (FL) and the reference layer (RL). 



III. EXPERIMENTAL RESULTS 

In this section, we describe the characteristics of the spec- 
tra measured for easy axis and hard axis applied fields as well 
as their dependence on the pillar size. Note that there is no 
device possessing simultaneously on EA and HA all the prop- 
erties stated as typical of a particular pillar size. The prop- 
erties described in the text are therefore those observed on a 
majority of the EA spectra and a majority of the HA spectra, 
but not always for the same devices. For the basic identifica- 
tion of the observed spin wave modes we recall in Fig. QI a X 
(c) the mode dispersion for a nanopillar consisting of a SAF 
and an ideal free layer (i.e. the latter is assumed not to in- 
teract with the SAF) in macrospin approximation. EA and 
HA spectrum each contain two types of modes: the uniform 
FL modes and the uniform acoustic mode of the SAF. The 
acoustic SAF mode is thereby the lower of the two SAF eigen- 
modes and corresponds to oscillations, for which the in-plane 
components of the SAF layer magnetizations stay antiparallel, 
i.e. the dynamical magnetizations oscillate 180° out-of-phase. 
The high-frequency second eigenmode of the SAF, the optical 
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FIG. 1: Macrospin description of a nanopillar of size S: frequen- 
cies of the uniform modes versus external field along (a) easy axis 
and (c) hard axis, calculated in macrospin approximation using 
as material parameters literature bulk or thin film values (see Ap- 
pendix IB 2b : magnetization 2.2 T for CoFe, 1.9 T for annealed 
CoFeB, exchange bias J eb = 4.5 x 10~ 4 I/m 2 , interlayer exchange 
jmt _ _g x iq~ 4 j/ m 2 ; the shape anisotropy fields were calcu- 
lated using demagnetizing factors extracted from OOMMF simula- 
tions. Panels (b) and (d) show the corresponding calculated hystere- 
sis loops. In panel (a), filled symbols are used for ascending field 
(P— > AP^ SF) and open symbols for descending field (SF— s- AP— > 
P). 



mode (in-plane components of dynamical magnetizations in- 
phase), is not detected in the measured frequency range up to 
26 GHz, and will therefore not be mentioned further. 



4 




1500 
900 



105 x 190 nm 2 ^ 




i i i i I i i i i I i i i i f^i i i | . . i ! 




— 




EI 









I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 



-150 -100 -50 50 100 150 200 
magnetic field \i H (mT, easy axis) 




N 

o 12 
c 

Q) 

§■ 8 





1 _ " - * _ 




10 

""I - 


-." -- - - • 












900-' 



-150 -100 -50 50 100 150 
magnetic field |J H (mT, hard axis) 



FIG. 2: Power spectrum density (PSD, log scale) versus descending FIG. 3: Power spectrum density (PSD, log scale) versus (ascending) 

(SF— >-AP— >P) easy axis applied field for a device of (a) size S, (c) hard axis field for the same devices as in Fig. [2] (a) size S, (c) size M, 

size M, and (e) size L. Panels (b), (d), and (f): corresponding hys- and (e) size L. Panels (b), (d), (f): corresponding hysteresis loops, 
teresis loops. 



A. Eigenexcitations for easy axis applied field 

1. Size independent properties 

For all pillar sizes, EA spectra (Fig. [2]) contain two groups 
of modes. The first group consists of V-shaped modes 
basically symmetric about zero field for high positive and 



negative fields and history dependent in the hysteretic field 
region with a discontinuity at the coercive field. Since this is 
the typical behavior of FL modes (Fig. Q2 a X (b)), the modes in 
this group are labeled with F. The second mode group consists 
of modes having a minimum at or near the spin-flop field of 
the SAF. This being the characteristics of the acoustic modes 
of the SAF (Fig. 02 a), (b)), the modes in this group are labeled 
with A. Modes that cannot be assigned unambiguously to 
one of the groups are labeled with U (like unidentified). 
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Within each group the modes are numbered consecutively 
with increasing frequency. The lowest FL mode FE has 
been labeled differently, because it shows in several aspects 
a qualitatively different behavior than the next higher mode 
FO. As we shall see in the following subsections, there is 
evidence that it might belong to excitations localized near the 
layer edges (hence the second label E). Occasionally, second 
harmonics, such as the mode 2FE in Fig. [2a), are observed. 



There are three FL modes common to all pillar sizes: the 
modes FO and FE, visible in the P and low-field AP state, and 
- with a size invariant spacing of 5 — 6 GHz to FO - the mode 
F3, mainly visible at high fields. 

In the SF region at high positive fields, several SAF modes 
are observed, the intense lowest modes having either one 
single minimum (Fig. [2je)) or two minima at different fields 
(Fig- Ha), (c)). The minima are always positioned in the 
vicinity of significant changes in the slope of the resistance 
in the hysteresis loop. The occurrence of more than one 
such slope change clearly indicates that the magnetization 
of the reference layer does not undergo a single abrupt 
transition as associated with the SF, but that there are several 
domains with different transition fields. Independent of 
the number of observed minima at the SF, the SAF modes 
show complicated, irregular structures, indicating strong 
non-uniformities of the SAF layer magnetizations in this 
field region, and will therefore not be discussed in more detail. 



The modes after the SF, which are labeled with U, comprise 
the lowest FL modes and higher order SAF modes. They 
cannot be identified with certainty, because both mode types 
have similar frequency, slope, and intensity in this field 
region. Moreover, above the SF, the large angle between 
the FL magnetization and the RL magnetization boosts the 
experimental sensitivity to both FL and SAF modes (see 
section IV Bi causing in particular for sizes S and M abrupt 
changes in the mode intensity, such that FL modes suddenly 
appearing above the SF field may be misinterpreted as SAF 
modes. 

For some devices, a change of mode intensity is also observed 
at a high negative field value (cf. FO, F3 in Fig. Et a ) at 
-140 mT, or AO in Fig. Etc) at -105 mT), at which the 
resistance changes by 10-20 Ohm (not visible on the scale 
of Figs. EtbX(d))- The reason for this is a change in the 
micromagnetic configuration, very likely of the RL (see 
section lVBT i. 



Finally, in the AP state, gaps with pillar size dependent 
opening from 2 to 0.5 GHz are observed in the modes FE and, 
if visible, F0 (see e.g. Figs. Et e ) an d ( a ))- We will see that 
this is a consequence of the mutual dipolar coupling between 
the FL and the SAF leading to anticrossing of FL and SAF 
modes. 



2. Size dependent properties 

In the following, we describe the size specific properties of 
the EA spectra. We will see that with increasing pillar size, 
the spectra undergo characteristic changes, some of which are 
caused directly by the increasing dimensions, while others 
are most likely consequences of an increasing non-uniformity 
of the magnetization. 

Size S 

In the EA spectra of the smallest pillars (Fig.Et a )) the modes 
FE, F0, and F3 are mostly the only observed FL modes. 
The opening of the anticrossing gap in F0 is typically 2 
GHz. F3 has for all devices a frequency of 22 — 23 GHz at 
+240 mT. FE and F0 show a slightly larger device-to-device 
variation: For some devices, the mode FE is visible up to very 
high negative fields (as in Fig. Et a )) with a roughly constant 
spacing to F0 of about 1.5 GHz, for others it is observed in 
the low-field region only where the resistance departs from 
its saturation value, i.e. where the (FL) magnetization shows 
signs of increasing non-uniformity; in this case, FE rapidly 
approaches F0 for increasing field and vanishes, once the 
resistance has reached its saturation value. This correlation of 
FE to the non-uniformity of the static magnetization suggests 
that FE might be an edge mode. Another observation in 
favor of this supposition is that the mode FE has a noticeably 
higher slope than F0 (due to its tendency to approach F0 
asymptotically from below), which would not be the case if 
both modes were volume modes. The average frequency of 
F0 is about 16 GHz at -190 mT and 6.5 GHz at zero field. 

SizeM 

The EA spectra of devices of pillar size M (Fig. Et c )) differ 
from those of the size-S devices in the following points: The 
mode FE has developed a minimum at low fields, which for 
some devices may almost reach zero frequency. This effect 
is likely to be caused by the increased non-uniformity of the 
magnetization in this field region. Occasionally, the mode 
F0, too, becomes deformed, though much less than FE. The 
relative intensity of F0 typically decreases, whereas that of 
FE increases. The observed evolution of FE with the pillar 
size, too, corroborates the assumption that this mode might 
be an edge mode. The frequencies of all modes decrease 
typically by 1 or 2 GHz at high fields. The gap opening in the 
modes FE and F0 has decreased to 1-1.5 GHz. At very high 
fields, additional modes, most likely belonging to the FL, 
appear just above or below the mode F3. 

Size L 

For the size-L devices (Fig. Ete)), the gap opening in the mode 
FE does not exceed 0.5-1 GHz, and the frequencies of FE and 
F3 have decreased by another GHz. The spacing between FE 
and F3 has not changed with respect to size M or S. Above F3 
three additional FL modes (F4 to F6) with a spacing of 1 GHz, 
and below F3 two extremely weak modes (Fl, F2) have ap- 
peared. As for size M, the mode FE is strongly deformed in 
the low field region and has still gained intensity with respect 
toFO. 
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B. Eigenexcitations for hard axis applied field 

1. Size independent properties 

In Fig.|3]are depicted the HA spectra of the same devices as 
in Fig. [2] The FL modes have a characteristic W-shape with 
two minima in the lowest modes at about ±70 mT. Though 
70 mT is for all pillar sizes substantially higher than the mea- 
sured anisotropy fields, the minima are often 20 interpreted as 
to correspond to the saturation of the free layer magnetization 
along the HA (cf. also Fig.[TJc),(d)). 

At or near zero field, the modes F0 and Fl show typically 
one, sometimes two sharp minima, which become deeper 
with increasing pillar size, and which are not present in the 
macrospin HA spectrum (Fig.QJc)). We therefore suspect that 
at least one of the modes observed at low fields is actually an 
edge mode FE. As a matter of fact, if at zero EA field the low- 
est mode is - as we think - an edge mode, then the lowest mode 
at zero HA field must be an edge mode, too, because zero EA 
and zero HA field are formally identical. The observed field 
dependence of the frequency of the mode F0 then implies that 
its character must be changing (continuously) from edge mode 
at low fields to volume mode at high fields. Such a progressive 
change of the mode character would be consistent with the ex- 
pected saturation process of the magnetization along the HA: 
for very low HA field, the magnetization in the central part of 
the layer is aligned along the EA due to the shape anisotropy, 
and only in narrow zones along the short edges of the rectan- 
gle the magnetization starts to align with the HA. For increas- 
ing field, these edge zones (domains) expand continuously to- 
wards the layer center, until for some field value the volume 
magnetization and finally the magnetization in the zones along 
the long edges of the layer saturate along the HA. The oscil- 
lations of this increasing part of the magnetization parallel to 
the field would obviously correspond to edge modes at low 
fields, and to volume modes at high fields; for intermediate 
field values, they would have a mixed character. 
Finally, for most devices, we also observe the almost horizon- 
tal lowest acoustic modes of the SAF, AO and Al (see e.g. Fig. 
|2a)), where AO has a frequency of typically 10 — 12 GHz at 
zero field. Note that in particular for size-S devices the SAF 
modes are strongly asymmetric w.r.t. zero field, even though 
the FL modes and the hysteresis loop are basically symmet- 
ric. As will be shown in section [Vll this asymmetry cannot 
be explained by a misalignment of the external field with the 
HA. 



2. Size dependent properties 

Size S 

The HA spectra of size-S devices are characterized by 
rounded saturation minima of the FL modes at ±70 mT (cf. 
Fig. [2a)). The frequency minima of the mode F0 are thereby 
not zero, but raised to values between 4 and 6 GHz, which is 
a consequence of the mutual dipolar coupling between the FL 
and the SAF, as we shall see in section [Vm Typically 5 to 7 
FL modes, F0 to F6 in Fig. [3 a), are observed, where the first 



two modes have frequencies of 12 and 14 GHz, respectively, 
at ±190 mT. The lowest mode F0 has much higher intensity 
than the other modes, and its frequency at high fields is the 
same for all measured devices of size S. The frequencies of 
the higher modes slightly vary from device to device. (In 
Fig- Sa) the modes F2 and Al seem to accidently coincide 
at positive fields; however, for other devices, F2 is clearly 
resolved.) The presence of Fl and F2 in the HA spectra for 
several devices of size S with rather high intensity is of great 
importance, because at least the mode F2 is not observed in 
the EA spectra, not even after the SF where the experimental 
sensitivity is comparable to that on the HA (see section [VBT l. 

SizeM 

In the spectra of pillar size M (Fig. 0c)) the saturation minima 
at ±70 mT in the mode F0 are much deeper than for size S, 
which is consistent with the lower dipolar coupling between 
the FL and the SAF concluded from the EA hysteresis loops 
and spectra (see also section [VTlb . The minima at zero-field 
are considerably sharper than for size S indicating increasing 
importance of edge domain effects. The acoustic SAF modes 
have with about 10 GHz at zero field for AO a slightly lower 
frequency than for size S, which is either due to a smaller 
interlayer exchange or the increasing non-uniformity of the 
SAF layer magnetizations (resulting e.g. on the EA in the 
observed stepwise switching in the SF region). Finally, the 
overall mode spacing has noticeably decreased compared to 
size S, as should be expected. 

Size L 

In the spectra of pillar size L (Fig. 0e)) the minima in the 
modes F0 and Fl at ±70 mT reach, as for size M, markedly 
lower frequencies than for size S. Both F0 and Fl are strongly 
deformed in the vicinity of their minima and may even cross 
each other. The minimum at zero-field has still become 
slightly deeper, the impact of edge domains now being domi- 
nant. Contrary to size S, the shape of the modes in the low and 
medium field region is strongly device dependent and sensi- 
tive to small changes of the field direction. 



IV. OUTCOMES AND LIMITS OF THE MACROSPIN 
MODEL 

Before making a detailed and rigorous analysis of the field 
dependence of the modes frequencies in the next section, 
we start by attempting to model the free layer modes F0 
using conventional Kittel fits. The aim is twofold: motivate 
the need for a more elaborate analysis by showing quantita- 
tive and qualitative limits of the macrospin approximation, 
and obtain approximate starting values for the magnetizations. 

Approximating the free layer as an isolated rectangular 
platelet with only shape anisotropy, its ferromagnetic reso- 
nance frequency is described by the well-known Kittel law, 
which for EA applied field (x-direction) reads 

uj 2 = j 2 [H appl + H k ][H appl + {N z - N X )M S ], 
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and for HA field (y-direction) 

uo 2 = ^[H appl - H k ][H appl + (N z - N y )Ms], 

where Ms is the saturation magnetization of the free layer 
and Hk = (N y — N x )Ms the shape anisotropy field. 
Applying Kittel fits to the modes FO in the high field regions 
of the spectra in Figs. |2][3] allows us to extract Ms and 
Hk for the different pillar sizes, independently for EA and 
HA. Using the demagnetizing factors N x , N v , and N z of 
Ref. [l5l we obtain from the modes FO in the EA spectra the 
following values for Ms and Hk- for size S HqMs = 1.14 T 
(for most devices of size S: 1.3 T) and /io-fffc = 37 mT, for 
size M hqMs = 1.04 T and /.loHk = 35 mT, and for size 
L hqMs — 0.91 T and noHk — 30 mT. In comparison, 
for the modes FE, larger magnetizations (1.3 to 1.1 T), 
but much smaller anisotropy fields (less than 12 mT) are 
obtained. Similarly, the modes FO on the HA yield for size S 
HqMs = 1.41 T and fiaHk = 55 mT (universal for size S), 
for size M /i Ms = 1.40 T and ^oH k — 76 mT, and for size 
L Li M s = 1.41 T and ^ H k = 78 mT. 

The minimum requirement for these values to be reason- 
able approximations is that the magnetizations and shape 
anisotropy fields extracted from EA and HA spectrum of the 
same device are roughly equal. However, as can be seen, both 
Ms and Hk are considerably larger on the HA, the discrep- 
ancies becoming larger with increasing pillar size. In addi- 
tion, on the HA - and, if the mode FE is used, also on the 
EA - the anisotropy fields are neither consistent with the ex- 
tracted magnetization nor with the anisotropies found by as- 
troid measurements^ (cf. section ITU). Therefore, treating the 
free layer and the SAF as uncoupled systems consisting of 
uniformly magnetized layers is obviously insufficient to de- 
scribe the eigenexcitations of nanopillars. The next section 
will be dedicated to a rigorous treatment of spin waves in a 
coupled three-layer system with lateral confinement. 



V. MODEL OF SPIN WAVE MODES IN NANOPILLARS 
A. Dipolar-exchange spin waves with quantized wavevectors 

1. Eigenexcitations of coupled three-layer system 

An MTJ nanopillar consists basically of three confined 
magnetic layers: the free layer, which will be labeled with 
the index "F", and below the two SAF layers - the reference 
(top) layer and the pinned (bottom) layer - labeled with in- 
dices "1" and "2", respectively. The magnetization dynamics 
in each layer I 6 {F, 1,2} of this coupled three-layer system 
is governed by the Landau-Lifshitz equation. 
For small amplitude precessions, the magnetization Mi (r , t) 
can be decomposed in zeroth order approximation into a time 
independent uniform (U) equilibrium component M ; (satu- 
ration magnetization Mi) and a small perpendicular dynami- 
cal part SMj 1 (r, t). Static non-uniformities of the equilibrium 



magnetization will be discussed in section [VB] The dynam- 
ical component SM^ (r, t) can be approximated as a sum of 
plane spin waves, 

SM^(f, t) = ^^2 ' (ky^-^n* , (i) 

k 

where the wavevectors k = (k x ,k y ,0) of the partial waves 
are quantized due to the spatial confinement of the layers. 
The out-of-plane component k z is zero for all modes in the 
experimental scope due to the very small layer thicknesses of 
2 — 3 nm. The quantization of the in-plane components k x , k y 
will be discussed in detail later on. The frequencies oj^ of 
the partial waves are the eigenfrequencies of the three-layer 
system. 

In the effective fields acting on the magnetizations the fol- 
lowing interactions have been taken into account: the applied 
field H appl , the exchange bias field acting on the bottom layer 
of the SAF (coupling constant J eb ), the interlayer exchange 
coupling of the SAF layers (coupling constant J mt ), and 
the (intralayer) exchange interaction in each layer (exchange 
stiffness constant Ai), as well as the demagnetizing fields and 
mutual dipolar coupling of the layers. For the demagnetizing 
fields we use the standard tensor expression for uniformly 
magnetized ellipsoidal bodies, where the diagonal compo- 
nents of the diagonal (self-)demagnetizing tensors N/ are the 
demagnetizing factors Nf, Nf ' , Nf of the rectangular layers. 
Although this approximation is expected to be satisfying 
for the static demagnetizing field, it is rather crude for 
the dynamical part, since the dynamical magnetization is 
non-uniform unless k = 0. 

The fields resulting from mutual dipolar coupling are given 
by analogous expressions where the (self-)demagnetizing 
tensors of trace 1 are replaced by the mutual demagnetizing 
tensors 21 N m / of trace (I, to G {F, 1, 2}, I ^ to). For the 
given pillar geometry, N m / is diagonal, too, as can easily 
be shown using the formulae for the tensor components in 
Ref. Hj]. The diagonal components will be referred to as the 
mutual dipolar coupling constants iV* j, NIL,, and N* nl . 
Note that there is no significant perpendicular surface 
anisotropy at the top and bottom surfaces of the layers in 
MTJs, as has been demonstrated in Refs. Il5ll22l Since the 
impact of the bias current on the experimental spectra has 
been found to be negligible, we do not include current-based 
interactions, such as spin-torque or the Oersted field. The 
latter does e.g. not exceed 1 mT for a current of 0.3 mA and 
an impact diameter of 100 nm. 

With these approximations the Landau-Lifshitz equations 
of the three pillar layers become a system of 3 x 3 = 9 coupled 
linear equations for the components of the dynamical magne- 
tizations SM^ (k). It can be solved as the eigenvalue problem 
of the 9x9 coefficient matrix F of the 9-component vector 
{5M%(k), 5M?(k),5M¥{k)) describing the dynamics of the 
three-layer system as a whole. The eigenvalues of F are the 
eigenexcitations uj^ of the three-layer system and can be cal- 
culated numerically as a function of the applied field, yielding 
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the expected spin wave spectra uj^(H appl ) of the nanopillar. 



2. Quantization of in-plane wavevector 

The in-plane components k x , k y of the wavevector are de- 
termined by the boundary conditions (BC) imposed on the 
dynamical magnetization (Q~|i at the lateral layer boundaries 
x = ±L x /2 and y = ±L y /2. For simplicity we will consider 
the x-component (along the long edge of the rectangle) as an 
example, where any of the following statements hold equally 
for the y-component with x and y permuted. 
For the x-component the BC read: 



= 0. (2) 



where £ x = x/L x . Eq. (O is a modified version of the 
effective BC derived by Guslienko et al. 9 for thin magnetic 
stripes. In difference to Ref. |9| we allow for different pinning 
parameters d+ and d x at opposite boundaries x = ±L x /2 
to account for potential asymmetries in the pinning expected 
from a real device. Moreover, instead of using the analytical 
expression (5) in Ref. HI to calculate the (dimensionless) 
pinning parameters, we will extract approximate values for 
d x from the experimental spectra (see section |Vl|i. 



Applying the BC (|2|i to the sinusoidal mode profile 

Jie lA:r = sin(k x x + tf> x ) sxn(k y y + (j> y ), 



(3) 



of the partial spin waves in SM^ (Q} yields for the wavevector 
component k x and the phase cj> x the quantization conditions 



=F k x L x cot I ±fc : 



L, 



] — Oj x 



(4) 



It is convenient to express k x L x in the argument of the cotan- 
gent as multiples of ir, thus defining the - in general noninte- 
gral - mode numbers 



k x L x 



(5) 



of the quantized spin wave modes (n x , n y ). 
For symmetric pinning, d x = d x = d x , it follows from 
that the cotangent has to be antisymmetric, yielding 
4> x = 7r/2 or 4> x = 0, i.e. symmetric or antisymmetric 
wavefunctions Q. In the limiting case of totally unpinned 
BC, d x = 0, the mode numbers n x are integers, starting 
at 0, and the corresponding wavefunctions alter between 
symmetric and antisymmetric for successive mode numbers, 
starting with symmetric, such that there are always antinodes 
at both boundaries. 

For finite values d x > of the pinning, the mode numbers 
n x are no longer integers. Plotting n x versus d x by means of 
eqs. and (O shows that with increasing d x , the deviations 
An x of n x from the corresponding integral values n x of 
the unpinned case increase continuously from An x = 



for d x = (unpinned) to An x = 1 for d x — oo (totally 
pinned). Therefore, the mode numbers for total pinning, 
— n x + 1, are integers again. For a fixed intermediate 
value d x , the deviation An x of the mode number n x from the 
corresponding integral mode number n x is found to rapidly 
decrease with increasing n x . For a given pinning, the mode 
numbers are therefore no independent variables: once one 
mode number (e.g. that of the lowest mode) has been fixed, 
all other mode numbers are fixed, too. 

In case of slightly asymmetric pinning, d+ ^ d~, the phase 
<j> x differs from the values (f> x ,a by a small phase shift A<j> x , 
such that the wavefunctions are no longer totally symmetric 
or antisymmetric. In this case, the mode numbers n x are 
necessarily non-integral. In the hypothetic case of totally 
asymmetric pinning, d+ = and d~ = oo (or vice versa), 
An x = 0.5 and A(j) x = A$" x %ax = tt/4. For arbitrary 
pinning, Acf> x is an unknown function of d x and n x . 

The pinning for a given in-plane direction of a magnetic el- 
ement depends on its dimensions and in addition on the inho- 
mogeneity of the internal field, 4 ' 9,23 Consequently, the mode 
numbers are expected to be larger for the a>direction than for 
the y-direction of the same pillar, and possibly different for 
easy and hard axis applied field. 



B. Expected experimental sensitivity 

1. Formulation of the problem 

As described in section|II] the experimental spin wave spec- 
tra are obtained by measuring the voltage noise of the pillar. 
To be more precise, we measure the average of the local volt- 
age noise over the pillar area. The local voltage noise is the 
product of the local current density and the local magneto- 
resistance (MR) noise generated by spin waves in the free 
layer (FL) and the reference layer (RL). In the ideal case of 
a homogeneous in-plane distribution of the current, the mea- 
sured voltage noise is proportional to the average of the lo- 
cal MR noise. For the sake of simplicity, we will derive the 
expected MR noise for excitations in the FL, where the anal- 
ogous expressions for the RL are obtained by permuting the 
indices "F" and "1". The consequences of inhomogeneities 
will be discussed later in this section. 

The MR noise signature of a partial spin wave with wavevec- 
tor k representing the FL mode (n x ,n y ) is in linear order 
given by the square of 
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SR F (k) = -— I Mtir) ■ 6M%(k) Ke lfc? dr (6) 
Opil Js p u 

where S p u is the pillar area, Mi (r ) the micromagnetic equi- 
librium magnetization of the RL, SMp(k) the amplitude vec- 
tor and 5fte lfer the spatial dependence (wavefunction) of the 
spin wave in the FL (cf. eq. (Q])). Decomposing Mi(r) as 
before into a uniform macrospin component , dominating 
in the central (volume) part of the layer, and the remaining r- 
dependent edge domain components 8M^(r), yields as final 
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(a) 



field 



easy axis 



below 2 nd SF 
P state 

AP state 
above SF 



H < HsF2 

highH 
low H 

low H 
highH 

H > H SF 



hard axis 



\H\>Q 
H « 



M?{6) 



Mi sin do 

M t A6 
Mi AO 

M t Ad 
M t A6 

Mi sin 9o 



Mi sin O 
Mi AO 



SM t E 



SM t E 



5M§ Sip 



SM? Sn 



8M? 



(b) 



integral W x (n x , 



mode 


for weak pinning 


for strong pinning 


number n x 


(An, < 1) 


(An, « 1) 





1 


2/tt 


odd 


A(j) x /n x ■ 2/tt 


A<j) x (l - An x )/n x 


even 


An x /n x 


l/n x ■ 2/tt 



TABLE I: Dependence of the magneto-resistance noise @ on the 
static micromagnetic configuration of layer / £ {F, 1} and the mode 
character, (a) leading order contributions of volume magnetization, 
M; (0), and edge domain contributions, 5Mf(n x , n y , <f> x , cj> y , 0), 
versus easy and hard axis applied field in different field regions. 
Si m is the Kronecker symbol, (b) leading order terms of integral 
W x (n x , <j) x ) of the wavefunction versus the mode number n x — 
n x + An, in the regime of weak and strong pinning. W v is given 
by analogous expressions. 

expression for the MR variation © 

SR F (n x ,n y ,4> x ,4> v ,9) 
= 5M$(k) [MY{6)W x {n x ,ct> x )Wy(n v ,<t> v ) 



x , (f>y,0) 



(7) 



The first term, MY W x W y , is the contribution of the uni- 
form volume magnetization to the MR noise (hence the su- 
perscript V). MY (9) = sin 9 Mi is the projection of M± onto 
SMp(k), where 9 denotes the angle between the macrospins 
of FL and RL. W x and W y are the integrals of the x- and y- 
dependent factors of the wavefunction (0, respectively. 
The second term, 5M^, is the contribution of the static edge 
domain magnetization components of the RL to the noise sig- 
nature of the (volume) FL mode (n x ,n y ) (not to be confused 
with edge modes). It is the spatial average of the projection of 
8MY{r) onto SMp(k) weighted by the wavefunction <j3j. 
Mathematical expressions for W x (W y ) and SM^ as well as 
details on the derivation of eq. (0 can be found in appendixlAl 
For the following discussion it is sufficient to consider the 
leading order terms of these quantities listed in Table U 
The terms MY and SM^, resulting from the equilibrium mag- 



netization, obviously depend on the static micromagnetic con- 
figuration of the pillar layers. The quantities W x , W y , and 

again SMf, involving the wavefunction, depend on the sym- 
metry properties of the spin wave mode. In the following, we 
will show which modes are expected to be observed in the ex- 
perimental spectra under which conditions by analyzing the 
micromagnetic configuration and mode character dependence 
of the experimental sensitivity. 



2. Micromagnetic configuration dependence of sensitivity 

The micromagnetic configuration of the pillar is sensitive 
to direction and strength of the external field, which is why 
in the following we distinguish between easy and hard axis 
applied field and identify field regions of distinct values M ; v 

and SMf (I G {F, 1}). 

The largest contributions to the MR noise and hence highest 
sensitivity are expected in the field regions of maximum 
volume magnetization contributions MY(9), i.e. where 
sin 9 ;§> 0. As can be seen from Table |Ja), a large zeroth 
order term sin 9q is found for non-zero HA fields, above the 
SF at positive EA fields and below the SF at negative EA 
fields. In the AP and P state, sin#o = 0, such that the leading 
order terms are of first order in A6>. 9q thereby denotes the 
angle between the FL and RL macrospins in case of an ideal 
pillar, and A9 a small deviation of 9 from 9q caused by a 
misalignment of the exchange bias field or the external field 
with the symmetry axes of the layers. 

The presence of edge domains on the EA, as marked in 
Table UJa), can be deduced from the hysteresis loops, as is 
explained in detail in appendix I A II Whether these edge 
domains give non-zero contributions SM^ to the MR noise 
depends on the symmetry properties of both the micromag- 
netic state and the wavefunctions (see appendix IA2I for 
details). 

In summary, high sensitivity to both FL and SAF modes can 
be expected on the HA at any finite field value and on the EA 
above the SF at positive fields. Weak higher modes will there- 
fore be visible, if at all, in these field regions (cf. Figs 12131 ). 
The sensitivity below the SF at negative EA fields is also en- 
hanced, though substantially less than for the other two cases. 
The change in intensity can be nicely seen in Fig.|2ta) at -140 
mT, and in Fig.[2|c) at -105 mT. 

For EA fields between the two SF fields, i.e. in AP and P 
state, modes become visible only through the misalignment 
A9 of the macrospins or through edge domains of appropriate 
symmetry. From the latter, slightly increased sensitivity is ex- 
pected for FL modes in the AP state at high positive fields just 
below the SF, and for both FL and SAF modes at low fields in 
P and AP state. The presence of edge domains may thereby 
entail the appearance of the corresponding edge modes in the 
spectra. 
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3. Mode character dependence of sensitivity 



The quantities W x , W y , and 5M® depend on the symme- 
try properties of the wavefunction of the mode (n x ,n y ). As 
the edge domain contributions become effective mainly in 
the low-field region, in which our model is in any case not 
expected to be accurate, we consider only the integrals W x 
and W y belonging to the volume magnetization contribution 
MY , dominating at high fields. 

In Table Q}b) the leading order terms of W x are listed as a 
function of the mode number n x = n x + An x in the regime 
of weak and strong pinning. n x , An x , and A<p x are thereby 
defined as in section IV Al The y-dependent factor W y is 
given by analogous expressions. 

For zero pinning (An x = A<f> x = 0), the fundamental 
mode n x = is the only visible mode. In the presence 
of pinning, the higher modes n x > 1 begin to appear: for 
symmetric pinning (An x > 0, A(f> x — 0) only those with 
symmetric wavefunctions (even n x ), in case of asymmetric 
pinning (An x > 0, A(f> x ^ 0) also those with antisymmetric 
wavefunctions (odd n x ). 

For weak pinning (An x <C 1), W x is for all higher modes 
n x > 1 of first order in a small quantity: in An x for even 
n x , in A(f> x for odd n x . In appendix lAl we show that the 
expected intensities of the higher modes (n°,n°) = (1,0), 
(0,1), and (2,0) are about two orders of magnitude lower than 
that of the fundamental mode (n°,n°) = (0, 0), whereas the 
mode (n°, n°) = (1, 1) is expected to have a four orders of 
magnitude lower intensity than (0,0). 

For strong pinning (An x w 1), the natural reference mode 
numbers are the mode numbers = n x + 1 of total pinning: 
even (odd) n x in the table correspond to odd (even) n^°. W x 
for symmetric wavefunctions (even n x , odd n^°) has now be- 
come a zeroth order quantity like for the lowest mode n x = 
(n^° = 1), whereas for antisymmetric wavefunctions (odd 
n x , even n^°) it has become second order in (1 — An x ) and 
A(f) x . Consequently, the lowest higher order modes close to 
riy* 3 ) = (3,1) and (1,3) will have intensities comparable 
to that of the fundamental mode (n^°, n^°) = (1,1). Modes 
with an even mode number n^y are expected to be at least 
four orders of magnitude weaker than (1,1). 

The above results on the expected relative mode intensity 
have been obtained under the assumption of homogeneous 
current density and homogeneous saturation magnetizations. 
Under these conditions the voltage noise is proportional to 
the MR noise, and - at high fields where the edge domain 
contributions are negligible - the MR noise is proportional to 
(W x W y ) 2 . As we have seen, in this case higher modes be- 
come visible if the integrals W x and W y are non-zero, that 
is for non-zero asymmetric pinning. However, even in the ab- 
sence of pinning, the measured voltage noise can be non-zero, 
namely if the saturation magnetization or the current distri- 
bution are inhomogeneous, because then the average (|6]i over 
the pillar area becomes an integral of a generally unharmonic 
- and for asymmetric inhomogeneities also asymmetric - func- 



parameter 




extracted value 


saturation 




1.27 T 




magnetization 


/-IQ1V12 


1.4 T 




exchange 


A 


18.0 x 10-^J/m 






14.0 x 10 _12 J/m 


exchange bias 


jeb 


4.5 x 10~ 4 J/m 2 


interlayer exch. 


lint 
J 


-3.9 x 10" 4 J/m 2 


lateral pillar 


L x 


100 ran 




dimensions 


Ly 


60 nm 




demagnetizing 
factors 


(N F ,N F ,N Z F ) 
(N^ 2 ,Nl 2 ,Nl 2 ) 


(0.035,0.065,0.9) 
(0.027,0.049,0.924) 


dipolar coupling 
constants 


(N F1 ,N F1 ,N F1 ) 
(N F2 ,N F2 ,N F2 ) 
(N? 2 ,N? 2 ,N? 2 ) 


(0.01,0.018, -0.028) 
(0.005, 0.009, -0.014) 
(0.007,0.012,-0.019) 


mode numbers 


(n x ,n y ) 


fOO, aOO 


(0.4,0.0) 


easy axis 




flO, alO 
fOl, aOl 
f20, a20 
fll, all 


(1.13,0.0) 
(0.4,1.0) 
(2.05,0.0) 
(1.13, 1.0) 


mode numbers 


(n x ,n y ) 


fOO, aOO 


(0.2,0.0) 


hard axis 




flO, alO 
fOl, aOl 
f20, a20 
fll, all 
f21, a21 


(1.05,0.0) 

(0.2,1.0) 

(2.02,0.0) 

(1.05,1.0) 

(2.02,1.0) 



TABLE II: Material and geometry parameters used to calculate the 
spectra and hysteresis loops in Fig. [4] Error bars for the parameters 
are given in the text. The labels fmn and amn associated to the 
mode numbers are those in Fig. [4] 



tion. For real devices, we may therefore expect finite sensitiv- 
ity to most of the higher modes. 



VI. EXTRACTION OF MATERIAL PARAMETERS 

In section IV Al we have derived the mode frequencies 
W£ as a function of the material parameters Mi, Ai, J eb , 
J mt , the geometry parameters (L x , L y ), (Nf, Nf , Nf), and 

C^mi' ^mV ^ml)' as we ^ as trie mo( ^ e numbers (n x ,n y ). 
In this section, we will finally extract these parameters 
from the experimental spectra. Since the model is based on 
the assumption of uniform equilibrium magnetizations, its 
application can be expected to be reasonable only for the 
smallest pillar size S, for which the non-uniformities of the 
magnetization had been found to be minimum. We assume 
that in the EA spectra, the modes F0 and F3, and in the HA 
spectra, the modes F0 to F5 and AO, Al are, at sufficiently 
high fields, volume modes describable by the model. 

As a matter of fact, not all of the above quantities are free 
input parameters to the model. In appendix|B]we show that on 
the basis of the measured layer dimensions, basic OOMMF 
simulations, previous works published in the literature, and 
a couple of reasonable assumptions, the number of free 
parameters can be reduced to the following quantities: the 
three mutual dipolar coupling constants Np lt Np 2 , Nf 2 \ 
magnetization Mp and exchange stiffness constant Ap of the 
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FIG. 4: Calculated mode frequencies versus magnetic field along 
easy axis (a) and hard axis (c) for a pillar of size S. Panels (b) and (d) 
show the corresponding calculated hysteresis loops. The parameters 
used to calculate the spectra are given in Tab.HH The modes f20, f 1 1 
have practically identical frequency, which is why only one mode is 
displayed. In panel (a), filled symbols are used for ascending field 
(P— > AP— ► SF) and open symbols for descending field. 



free layer; interlayer exchange J mt and exchange bias J eb of 
the SAF; the mode numbers of the lowest FL and SAF mode. 
In appendix|B] we extract minimum and maximum values for 
the remaining parameters by adjusting the calculated modes 
and hysteresis loops to the corresponding experimental 
data (Figs. I2l3r a).(b)). In particular, we show that from the 



experimental constraints it follows that the pinning of the 
magnetization at the boundaries must be weak. 
Best overall agreement of experiment and theory is obtained 
for the parameters in Tab. [TT] The calculated spectra and 
hysteresis loops are shown in Fig. [4] In the following, we will 
point out similarities as well as differences between experi- 
ment and theory, and will discuss the "technical" reliability 
of the parameter values. The physical consequences will be 
discussed in the next section. 

The HA spectra (Figs. Oa)|4fc)) show quantitative agree- 
ment for the FL modes F0 to F3 at medium and high fields, 
and for the lowest SAF mode AO at high negative field 
only. However, there is no calculated mode corresponding 
to the mode F4; the frequency of the mode til, which is 
the next higher mode after f20/fll, is much too high for 
F4. A possible reason for this discrepancy might be our for 
higher modes rather crude approximation of the dynamical 
demagnetizing field (see below). The systematic asymmetry 
of the SAF modes, resulting in a much lower frequency of AO 
at positive field, cannot be accounted for by a tilting of the 
field w.r.t. to the HA, as this would affect both FL and SAF 
modes as well as the hysteresis loop. It might rather be caused 
by a misalignment of the exchange bias field, determining 
the magnetic symmetry axis of the SAF, with the geometrical 
symmetry axes of the rectangle, coinciding with the magnetic 
symmetry axes of the free layer. 

In the zero-field region, the differences between calculated 
and experimental HA spectra become substantial, as the 
measured modes develop pronounced minima, whereas the 
model predicts a local frequency maximum. 
On the EA (Figs. |2j4fa)), the calculated modes fOO and 
f20/f 1 1 fit the average experimental modes F0 and F3 (see 
section IIII Ai rather well. In particular, the gap opening in 
the mode F0 is reproduced in the theoretical spectrum. The 
modes fOl (F2), and - if none of the unidentified modes U 
is Fl- also flO, are not observed in the experimental EA 
spectra. Their absence is likely to be due to either an overall 
lack of intensity, first noticeable for the weakest modes, 
or a lower sensitivity to these particular modes on the EA, 
although our considerations in section lVBl vield no satisfying 
explanation for the different visibility for EA (above the 
spin-flop) and HA field under the made assumptions. For the 
SAF modes there can be only qualitative agreement due to the 
multiple SF in the experiment, which - like the low-frequency 
supposed edge modes FE - can of course not be described by 
a macrospin-based model. 

Finally, the calculated hysteresis loops in Fig.Htb), (d) are in 
qualitative agreement with the measured loops. 

Main error sources for any of the parameters are obviously 
the various assumptions in section [VAl and appendix[B] Addi- 
tional uncertainties come from insufficient experimental data, 
e.g. the ignorance of the actual SF field or device-to-device 
diversity of mode frequencies. In particular, our represen- 
tation of the dynamical demagnetizing field in the standard 
tensor expression for uniformly magnetized ellipsoidal bodies 
can be expected to be a reasonable approximation only for the 
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lowest mode, because there the dynamical magnetization is 
indeed almost uniform. As we extract the material parameters 
mainly from this mode, their values are only little affected by 
this approximation. The largest discrepancies are expected for 
the exchange stiffnesses, as they rely necessarily on the higher 
modes. 

The value of fioMp is found to be between 1.25 T and 1.3 T , 
Ap is expected to lie in the interval (18 ± 3) x 10~ 12 J/m. 
Under the assumptions in appendix [B] the same holds true 
for Mi and A\. The error of M 2 and A 2 is significantly 
larger than for the FL and RL, because of the additional de- 
pendence on the thin film value of the CoFe layer and the 
ignorance of the experimental exchange bias field. H0M2 
is expected to be contained in the interval (1.4 ± 0.1) T, 
and A 2 in (16 ± 4) x 10~ 12 J/m. J mt is estimated to be 
(-4.0 ± 0.4) x 10- 4 J/m 2 , and J eb (4.2 ± 0.7) x 10^ 4 J/m 2 . 
The mode numbers of F0 (fOO) have to be smaller than 
(0.6, 0.0), or (0.3, 0.2) for n y > 0, in order to ensure sat- 
isfactory agreement in frequency and a reasonable value for 
Ap. The agreement is better, if the mode numbers are chosen 
smaller on the HA than on the EA, and n x > n y . For the sake 
of simplicity, we have therefore set n y = 0. 
N F1 = 0.01 and N F2 = 0.005 are uniquely determined 
by the experimental constraints with a maximum deviation of 
±0.002. The deviations of ±10 nm of the lateral dimensions 
L x , L y from the mean values will change all geometry related 
parameters accordingly. 



VII. DISCUSSION 

In the previous sections, we have modeled the spin wave 
spectra of MTJ nanopillars as eigenexcitations of a coupled 
three-layer system with lateral confinement. In this section, 
we will see, which properties of the experimental spectra can 
be explained in the scope of this analytical model, and which 
cannot. First, we will discuss the material parameters of the 
pillar extracted from the high field regions of the spin wave 
spectra. Thereafter, the low-field anomaly of the spectra and 
its relevance for applications will be discussed. Finally, we 
will summarize the properties of the experimental spectra, 
which are beyond the approximations of our model, includ- 
ing the pillar size dependence. 



A. Material and geometry parameters 

In this subsection, we will discuss the physical relevance 
of the extracted parameter values of the pillar in Table UTI 

With 1.27 T the saturation magnetization of the 
Co6oFe2oB2o layers of the pillar is significantly reduced 
compared to the thin film value of (1.8 ± 0.1) T^ 4 ^ or the 
bulk value for the underlying Co 75 Fe 2 5 of (2.2 ± 0.1) 
A reduction of the magnetization in nanopillars has already 
been observed in previous studies on pillar devices^— 
Three scenarios are usually suggested to account for this 
phenomenon: process-induced damages, 15 ' 32 current-induced 



heatingriai or a nonlinear change of the frequency with high 
mode amplitude^ 5 - As we work with low bias current, current 
induced heating can be excluded in our case. Similarly, 
since spin-torque induced auto-oscillations in our samples 
occur typically for currents above 1.6 mA for size L, 19 high 
amplitude nonlinear effects as possible cause can be rejected, 
too. Therefore, some sort of process damage, such as ion 
implantation, diffusion, intrinsic chemical modifications or 
interface effects, must be at the origin of the magnetization 
reduction, whose further investigation exceeds the scope of 
this paper. 

Concerning the boundary conditions and exchange stiffness 
we have come to the following conclusions: Strong pinning 
can be ruled out in our pillars (see section IVTb ; reasonable 
agreement between calculated and experimental data is 
obtained under the assumption of weak pinning. The pinning 
parameter deduced from the extracted mode numbers is with 
d < 1 about 10 times smaller than the one calculated by 
means of eq. (5) in Ref. 9 (d « 10) when using the material 
parameters of Tab. [II] Any value of d substantially larger than 
3 is found to yield mode numbers for the lowest mode very 
close to 1, i.e. strong pinning. This discrepancy between our 
result and the predictions of Guslienko's analytical model 9 
are not understood, as the latter is expected to be valid in 
the regime of element thicknesses smaller than the exchange 
length as well. 

We emphasize that, just as the magnetization, the exchange 
stiffness of the free layer does not exceed 2/3 of the thin film 
value, independent of the boundary conditions. Therefore, 
the magnetic properties of the nanopillar can by no means 
be described by the values measured on the unprocessed thin 
films. 

The mutual dipolar coupling accounts for several features 
of the experimental spectra: In the HA spectra, the mutual 
dipolar coupling of the FL and the SAF raises the frequency 
minima of F0 by several GHz, pushes them to slightly higher 
fields, and lowers the slope of the modes, reducing their 
frequency at ±190 mT by about 1 GHz. It also causes the 
bell shape of the HA hysteresis loop, by forcing the pillar into 
the AP state at low fields, and smoothes out the sharp bends 
at the anisotropy fields, which are observed in the case of an 
uncoupled free layer. 

In the EA spectra, the gap openings in the mode F0 stem 
from the anticrossing of F0 with the acoustic SAF modes 
due to coupling-induced mode hybridization. Finally, the 
net dipolar coupling field created by the SAF layers and 
favoring the antiparallel configuration of the pillar causes a 
shift of the EA hysteresis loops to negative fields of 5 mT, 
which is approximately 50% of the observed total shift. The 
remaining 50% may be due to an unequal reduction of the 
coercive fields at positive and negative field, which occurs 
if the micromagnetic configuration causes the FL to switch 
more easily from the P state to the AP state, than from the AP 
state to the P state. Indeed, the FL magnetization is expected 
to be more non-uniform - and consequently easier to switch - 
in the low-field P state because of the mutual dipolar coupling 
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field pointing antiparallel to the magnetizations in the P state, 
but parallel in the AP state. 

Within the diagonal tensor approximation of the mutual 
dipolar coupling, the tensor components are found to be 
significantly smaller than the values predicted by the for- 
malism developed by Newell et al.— or by the simplified 
version using for the in-plane components of the mutual 
dipolar coupling tensor the corresponding components of 
the self-demagnetizing tensor, as is commonly practiced 
when modeling flip-flop switching in MRAM celki^i The 
coupling between e.g. the free layer and the reference layer of 
our pillars would be overestimated by Ref.[2l]by a factor of 2, 
and by Refs. [3611381 by a factor of 3. A possible explanation for 
this reduction of the interlayer dipolar coupling may be that 
the coupling field extracted from the experiment is actually an 
effective mutual dipolar coupling field comprising the dipolar 
coupling due to the charges at the lateral layer boundaries as 
well as some Neel-type coupling resulting from the correlated 
roughness of the three magnetic layers. This orange-peel 
coupling may partially compensate the antiparallel coupling 
due to the charges at the layer edges. Another possibility is a 
reduction of the dipolar coupling due to the non-uniformity of 
the micromagnetic magnetization at the layer edges, though 
this effect should be small at high fields. 

Finally, the extracted exchange bias energy and the inter- 
layer exchange coupling are consistent with the large body of 
dedicated literature (see e.g. the values in appendix lB 21 ). 



B. Low-field behavior and its relevance for applications 

In sections [II] and [TTTJ we have seen that at low fields both 
EA and HA spectra show for all three pillar sizes unmistak- 
able signs of non-uniform magnetizations: in the HA spectra, 
the FL modes possess at zero field, instead of the local max- 
ima predicted by the model, sharp minima, whose depth in- 
creases with increasing pillar size, indicating increasing non- 
uniformity of the magnetization. In fact, the modes FO, Fl 
are likely to change character from volume modes at high and 
medium fields to edge modes at low field. The EA spectra 
contain low-frequency supposed edge modes FE, which be- 
come progressively deformed around zero field for increasing 
pillar size, i.e. for increasing non-uniformity of the magneti- 
zation. 

The non-uniformities of the magnetizations are expected to 
influence the switching dynamics of the pillar. The first con- 
sequence is that they lower the coercive field (as discussed in 
the previous paragraph), thus enlarging its difference to the 
shape anisotropy field. This effect has indeed been found to 
be particularly strong for size L (see section[II]l. More impor- 
tantly, the fact that the lowest mode is not the uniform mode, 
but an edge mode, will affect the magnetization reversal path 
in current-induced switching, favoring non-uniform reversal 
paths, as has already been concluded indirectly from reversal 
speed experiments^ 



C. Spin wave phenomena beyond the analytical 
approximations 

Based on the assumption of macrospin equilibrium magne- 
tizations, our model is certain not to describe any effect re- 
sulting from non-uniformities of the magnetization. This is 
the case e.g. for the low-field behavior discussed in the pre- 
vious paragraph, or the occurrence of more than one spin-flop 
transition of the SAF at positive EA field. 
However, even for high fields and pillar size S where the 
model is expected to work reasonably well, there are qualita- 
tive discrepancies between calculated and experimental spec- 
tra in frequency or visibility of higher order modes. 
Similarly, the high-field evolution of the spin wave spectra 
with the pillar size is not consistent with the predictions by the 
model. Although the model allows to reproduce qualitatively 
the EA spectra for pillar size L under reasonable assumptions, 
it fails for the high-field HA spectra. 



VIII. CONCLUSIONS 



In this paper, we have studied the magnetic field depen- 
dence of the mode frequency of thermally excited spin waves 
in rectangular shaped MgO-MTJ nanopillars of different lat- 
eral sizes. The spin wave spectra (frequency versus easy and 
hard axis applied field) of individual devices were obtained us- 
ing spectrally resolved electrical noise power measurements. 
In all spectra, several independent quantized spin wave modes 
stemming from eigenexcitations in the free layer and the SAF 
layers of the MTJ have been observed. By diagonalizing the 
dynamical matrix of a system of three coupled, spatially con- 
fined magnetic layers, we have modeled the mode frequencies 
for the smallest pillar size, 60 x 100 nm 2 , obtaining quanti- 
tative agreement for a majority of modes at high and medium 
applied fields. Our ability to detect a particular spin wave 
mode depends on the static micromagnetic configuration of 
the layers as well as on the symmetry properties of the mode. 
With the help of these discrimination criteria, we could iden- 
tify the observed modes and extract the material parameters 
of the pillar (Tab. Hill. The magnetizations and exchange stiff- 
ness constants were found to be significantly reduced com- 
pared to the corresponding thin film values, whereas the inter- 
layer exchange coupling and the exchange bias are consistent 
with their thin film counterparts. The interlayer dipolar cou- 
pling between the different layers could be well described in 
terms of an effective mutual dipolar coupling. Moreover, we 
could infer that the pinning of the magnetizations at the lateral 
boundaries must be weak. 

Finally, at low fields and for larger pillar sizes, there is clear 
evidence for strong non-uniformities of the layer magneti- 
zations, leading to qualitative differences between calculated 
and measured spin wave frequencies. 
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Appendix A: Derivation of magneto-resistance noise signature 

In this appendix we derive eq. (0 and the leading order 
terms in Table Q] starting from eq. ©. 

The first step is to evaluate the dot product in §6§ of the 
dynamical magnetization 5Mp(k) of the FL with the equilib- 
rium magnetization Mi(r) of the RL. Since SMp(k) is per- 
pendicular to the macrospin component Mp 7 , the dot product 
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of 8Mp(k) with MY can be expressed in terms of the an- 
gle 9 between the macrospins Alp 1 and Mf 7 of the two layers. 
Similarly, decomposing 8M^(r) into a component SM^,, (r) 

parallel to the macrospin Mj 7 , and a component SMf ± (r) 
perpendicular to , allows to evaluate the dot product of 
8Mf(r) with SM^(k). With that the MR variation © be- 
comes 



SR F {k,9) = SM^ik) MY{k,9) + Mf(k,6) 



(Al) 



where 



MY{k,0) = sin^Mj W x (n x ^ x ) W y {n y ,<p y ), (A2) 



1 



W x (n x , <t> x ) W y {n v ,4>y) = — / Me lkr df, (A3) 



and 

Mf {%,&) 
1 

Spil 



cos 9 6Mf ± (r ) + sin 9 5M®„ (r) 



(A4) 



As the edge domain contributions ( 1A4I > are relevant only on 
the EA in P and AP state where 9 is basically or it, the term 
with sin 9 in dA4b is in all practical cases negligible, such that 
only the term SAlYj_(r) = SAlY y (x, y) remains. 
In the following, we derive the leading order terms of these 
quantities listed in TableU 



1. Micromagnetic configuration dependence 

Decomposing 9 into the angle 9q between the two 
macrospins for an ideal pillar, and a small deviation A9 
due to misalignments, sin 6* and cos9 in MY(k,9) and 
AlY(k, 9) can be expanded in A9 about 9q, where the leading 
order terms for the different field regions are summarized 
in Tab. QIa). The underlying values of 9o are as follows: 
For EA applied field, 9q — in the P state, 9q = ir in the 
AP state, and -C Oq < ^ above the SF at positive fields, 
and below the 2nd SF at negative fields. For HA field, 9o 
decreases continuously from ir at zero field to a value close 
to 7r/2 at the saturation field of the FL, and finally towards 
zero as the RL magnetization continues to tilt towards the HA. 

The presence of edge domains on the EA can be deduced 
from the hysteresis loops: deviations of the resistance from its 
saturation values in P and AP state indicate non-uniformities 
of the FL and/or the RL magnetization. E.g. in the P state, the 
resistance increases continuously when the (ascending) field 
approaches the switching field to the AP state, both magne- 
tizations being subject to an increasing effective antiparallel 
field consisting of the (self)demagnetizing field, the mutual 
dipolar coupling field, and the external field as soon as it be- 
comes positive. 



In the AP state, at low negative field just before the switch- 
ing to the P state, both the external field and the interlayer 
dipolar coupling field are parallel to the magnetization of the 
RL, thus partly suppressing the edge domains created by the 
(self)demagnetizing field; in contrast, edge domains in the FL 
are only suppressed by the interlayer dipolar coupling field, 
but enhanced by the external field. Indeed, the resistance in 
the low-field AP sate departs much less from AP remanence 
than it does from the P remanence in the low-field P state. 
Finally, in the high-field AP state, the resistance starts to 
decrease continuously already long before the SF transition 
due to increasing non-uniformities of the RL magnetization, 
which is pointing antiparallel to the high external field. 
On the HA, edge domain contributions are negligible com- 
pared to the zeroth order contributions of the volume magne- 
tization, except for zero field where the macrospins are an- 
tiparallel. 



2. Mode character dependence 

By means of eqs. (O and Q the two integrals W x and W y 
over the layer dimension in direction x and y, respectively, are 
easily evaluated as 

2 7T 

W x (n x ,(j) x ) = sin(n x -) sin <fi x , (A5) 

n x ir 2 

where W y (n y , cf> y ) is given by an analogous expression. 
Decomposing the mode numbers and the phase as in 
section IV A| n x = n x + An x and <fi x = 4> x + A(f> x 
where <fi x = n/2 + n x ■ tt/2 is the phase for symmetric 
pinning, W x (n x ,(f> x ) can be expanded in A<fi x ^ 1 and 
either An x ^ 1 (weak pinning) or (1 — An x ) -C 1 (strong 
pinning). The result as a function of n x is shown in TableUb). 



The edge domain contributions SMf to the MR noise in 
the P and AP state can be evaluated on the basis of symmetry 
considerations. In spite of a non-uniform equilibrium mag- 
netization, SMp is zero if the product of the wavefunction 
and the function describing the spatial dependence of the 
y-component of the edge domain magnetization under 
the integral ( IA4l i is either zero or antisymmetric in x- or 
y-coordinate. For strong pinning, the product of the two 
functions is zero (or negligibly small), because near the layer 
edges, where the edge domain magnetization is non-zero, the 
wavefunction has minimum amplitude due to the pinning. 
Significant contributions from edge domains can be expected 
for weak pinning only. In this case, the integral (IA4I) will 
vanish for certain modes if the magnetization for a given 
micromagnetic state is invariant under reflection or rotation or 
a combination of both. The flower-state e.g. is invariant under 
reflections about x- and y-axis, i.e. the y-component of the 
magnetization is antisymmetric in both x- and y-coordinate. 
SMp is therefore non-zero only for modes with two odd 
mode numbers. Similarly, it can be shown that for the S-state, 
6 Alp is non-zero for modes, whose mode numbers are either 
both odd or both even; the C-state renders modes with odd n x 
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visible. 

Finally, we derive the expected relative intensity of the 
modes. An x for weak pinning or (1 — An x ) for strong pin- 
ning are of the order 0.1. We may assume that for small 
asymmetries of the pinning, Acf> x is at most of the same or- 
der of magnitude as An x (or (1 — An x )). Therefore, for weak 
pinning, W x (n x > 1) « 0.1 /n x is one order of magnitude 
smaller than W x (n x = 0) as 1, and we expect to observe 
in addition to the quasi-uniform mode close to (0,0) higher 
modes with mode numbers close to (1,0), (0,1), (2,0), (0,2), 
and possibly (3,0). Their intensities, being proportional to 
(W X W V ) 2 , scale with factors quadratical in A<j> Xiy or An XtV , 
and are therefore two orders of magnitude lower than that of 
(0,0). The intensities of all other modes, such as (1,1), are of 
forth order in A<f> x , y and An x , y , or strongly reduced by the 
factor 1 / '(n x n y ) 2 , and therefore most likely too weak to be 
detected. 

For strong pinning and even n x , W x (n x > 1) ~ l/n x is 
of the same order of magnitude as W x (n x — 0) ~ 1. In 
contrast, for odd n x , W x (n x > 1) w 0.01 /n x is two orders 
of magnitude smaller than W x (n x = 0). Consequently, the 
higher modes close to (3,1), (1,3), and (5,1) will have intensi- 
ties comparable to that of the fundamental mode close to (1,1), 
or one order of magnitude lower due to the factor 1 / (n x n y ) 2 . 



Appendix B: Details on extraction of model parameters 

In this annex, we present the arguments used to extract the 
material parameters Mi, Ai, J eb , J mt , the geometry param- 
eters (Nf, Nf , Nf), and (N^, N^, N^), as well as the 
mode numbers (n x , n y ) from the experimental spectra. 



constants A l : A 1 = A F and Af lar /A f 2 ilm = Af[ ar /A f *l 



2. Literature values 

In this paragraph, we list as an orientation literature values 
for the material parameters. 

As thin film exchange stiffness constants we use the values 
of the 40 nm CoFeB and CoFe films in Ref. H A f F u ™ = 

28.4 x 10- 12 J/mand A f 2 ilm = 27.5 x 10~ 12 J/m. 
The magnetizations of CoFe and annealed CoFeB depend 
on the percentage of Fe in Co: The bulk value for both 
Co7oFe3o and Co7sFe25 (corresponding to Co6oFe2oB2o) is 
(2.2 ± 0.1) T^i2^. The (thin film) free layer magnetiza- 
tion of our MTJ stack has been measured to be fXvM F %lm = 
(1.8 ± 0.1) TJ£ The thin film value for the CoFe layer is ex- 
pected to be in the interval /x M 2 /i/m = (2.0 ± 0.2) T. 
The exchange bias field in a CogoFeio (5 nm)/ PtMn (20 nm) 
system has been measured to be /io-ff eb ~ 67 mT^^Z, which 
corresponds to an exchange bias energy of J eb — 4.5 x 
10 -4 J/m 2 , using 2.0 T as saturation magnetization of the 
CoFe layer. For the interlayer exchange energy a maximum 
value of J mt = —5 x 10~ 4 J/m 2 has been reported. 28 
Ref. Hj] allows to calculate the dipolar coupling constant for 
two rectangular layers of equal thicknesses. As in our pil- 
lars the FL has a different thickness than the two SAF layers, 
only Nf 2 may be calculated directly, yielding N'f 2 = 0.016. 
The dipolar coupling constants N F1 and N F2 involving the 
FL can only be estimated as the mean value of the constants 
calculated for two 3 nm thick layers and for two 2 nm thick 
layers, from which we obtain N F1 « 0.018 and N F2 « 0.013 
(maximum deviation ±0.003). 



1. Reduction of number of free parameters 

Given the (approximate) layer dimensions L x ,L yi L z , the 
demagnetizing factors Nf, Nf, Nf can be calculated using 
e.g. OOMMF simulations, where we find Nf w 1 — (Nf + 
Nf) and Nf /Nf f» L x /L y as should be expected. Using 
the formulae in Ref. |2l| it can be shown that the dipolar cou- 
pling constants obey similar relations, N]L,/N*, = L x /L y 
and Nf nl = —(N^ + N^), and for symmetry reasons, 
N m ; = N/ m . The remaining components N F1 , N F2 , and 
Nf 2 are kept as free parameters to be extracted from the exper- 
iment, although they can be calculated by means of Ref. ml 
On the basis of previous measurements on MTJ stacks, the 
number of free parameters can be further reduced: In Ref. 22 
it has been shown that the magnetization of the CoFeB free 
layer does not depend on the layer thickness in the range from 
2 to 3 nm. We may therefore assume that the FL and the RL - 
being of the same material, but having different thicknesses - 
have equal magnetizations, Mi = Mp. Moreover, we expect 
the layer magnetizations in the pillar to be reduced for all lay- 
ers by the same (relative) amount w.r.t. the thin film saturation 
magnetizations: Mf llar /M{ ilm = M p F J ar '/M f F u ™ ■ Analo- 
gous relations are expected to hold for the exchange stiffness 



3. Regression method 

Mp is determined by the modes F0 on the EA with a 
weak dependence on the chosen mode numbers of F0 (see 
discussion below). A minimum value for Mp of 1.25 T fol- 
lows from the measured room-temperature anisotropy field, 
which must be smaller than the calculated (zero-temperature) 
anisotropy field. M\ and M 2 cannot be extracted directly, 
but depend entirely on the above assumptions. Once M\ and 
M 2 have been fixed, J mt and J eb can be estimated from the 
spin-flop field and the mode AO on the HA. N F1 and N F2 
follow from the gap opening in the mode F0 on the EA and 
the shift of the EA hysteresis loop to negative fields: We have 
used the above calcuted values for N F1 , N F2 , and A^f 2 as 
starting values, which we have adapted to the experimental 
data by rescaling, assuming that the deviation of the mutual 
dipolar fields, e.g. due to micromagnetics, is similar for all 
pillar layers. It turns out that N F1 as the largest component 
can be maximum 0.01, because otherwise the gap opening 
exceeds the observed 2 GHz (maximum value of 2.5 GHz for 
the calculated coefficients). On the other hand, the difference 
of N F1 and iVJ^ must be at least 0.005 to ensure a shift 
of the hysteresis loop of minimum 5 mT. N F1 — 0.01 and 
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Np 2 — 0.005 are therefore uniquely determined (maximum 
deviation ±0.002). Since for the calculated constants Nf 2 is 
between N F1 and N F2 , we set N'f 2 = 0.007. 

The exchange stiffness constant Ap and the mode num- 
bers (n x , n y ) cannot be extracted separately, since they enter 
the effective field (and consequently the frequencies) only as 
a product. The BC in the pillar being unknown, the lowest 
modes can have any mode numbers between (0,0) (unpinned 
BC) and (1,1) (totally pinned BC), where n x can be larger 
than n y (cf. section [VA] i. 

In order to adjust F0 on the EA with the mode (1,1) in the 
limit of strong pinning, we would need Ap » 1 /20-A F tlm and 
HoMp w LIT; this value for Mp is significantly smaller than 
the allowed minimum, and the reduction of Ap w.r.t. its thin 
film value is unreasonably large given thatM F w 2/?>-M f F am . 
In addition, a discrepancy of more than 1.5 GHz between cal- 
culated and measured mode F0 is observed on the HA even 
at high fields. Similarly, fitting F0 with the mode (0.5,0.5), 
which might be considered as the border between strong and 



weak pinning for n x — n y , or the mode (1,0), for maxi- 
mum difference n x — n y , still requires Ap < 1/5 • A f F Um 
and fioMp w 1.2 T. Consequently, strong pinning can be ex- 
cluded in our pillars; the mode numbers of F0 must be well 
below (0.4,0.4) or (0.8,0). The pinning is weak. This is also 
corroborated by the fact that in particular on the HA the low- 
est mode F0 has much higher intensity than the higher modes 
Fl to F5, which is a characteristics of weakly pinned systems 
(cf. section fVBT i. 

If we assume totally unpinned BC - fitting F0 with (0,0) and 
the higher modes with (1,0), (0,1) etc. -weget^oMp = 1.3T 
and Ap ss 2/3 • A p m , i.e. approximately Ap cc Mp. Tonar- 
row the mode numbers down within these borders, we assume 
that indeed n x > n y , which finally confines (n x ,n y ) on the 
HA to (n x ,n y ) < (0.2,0.1), and on the EA to (n x ,n y ) < 
(0.4, 0.2). To fit F0 on EA and HA simultaneously, requires 
that the mode numbers on the HA are smaller than on the EA, 
as had already been suggested in section [VAl 



